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ABSTRACT 

We have discovered 11 ultra-faint (r < 22.1) low surface brightness (LSB, central surface brightness 
23 < /x,- < 26) dwarf galaxy candidates in one deep Virgo field of just 576 arcmin 2 obtained by the 
Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT). Their association with the 
Virgo cluster is supported by their distinct position in the central surface brightness - total magnitude 
plane with respect to the background galaxies of similar total magnitude. They have typical absolute 
magnitudes and scale sizes, if at the distance of Virgo, in the range —13 < M r < —9 and 250 < r s < 850 
pc, respectively. Their colors are consistent with a gradually declining star formation history with a 
specific star formation rate of the order of 10 -11 yr” 1 , i.e. 10 times lower than that of main sequence 
star forming galaxies. They are older than the cluster formation age and appear regular in morphology. 
They represent the faintest extremes of the population of low luminosity LSB dwarfs that has been 
recently detected in wider surveys of the Virgo cluster. Thanks to the depth of our observations we 
are able to extend the Virgo luminosity function down to M r ~ —9.3 (corresponding to total masses 
M ~ 10 7 M 0 ), finding an average faint-end slope a ~ —1.4. This relatively steep slope puts interesting 
constraints on the nature of the Dark Matter and in particular on warm Dark Matter (WDM) often 
invoked to solve the overprediction of the dwarf number density by the standard CDM scenario. We 
derive a lower limit on the WDM particle mass > 1.5 keV. 

Subject headings: cosmology: observations — galaxies: clusters: individual (Virgo) 


1. INTRODUCTION 


The predicted over-abundance of galaxy satellites is a 
well known challenge to the standard Cold Dark Mat¬ 
ter (CDM) scenario, and it appears to be related to 
the large small-scale power of the CD M power spectr um 
of density perturbations (see e.g., [Klypin et al. I (|1999l )). 
The inclusion of baryonic processes in the models can re¬ 
duce this discrepancy (e.g., by suppressing star formation 
through supernovae feedback or photoionizing UV back¬ 
ground, [Governat^etjir] (l2007f) j. although it is still diffi¬ 
cult to maintain a simultaneous matching with other prop¬ 
erties of satellite galaxies, like colors and star-formation 
histories. An interesting alternative is to assume that 
Dark Matter is composed of p articles with mass wy ~ 1 
keV (warm Dark Matter, e.g.. [Sommer-Larsen fc Dolgov 1 
(12001ID . The larger corresponding free-streaming length 
results in a strong supp ression of the powe r spectrum at 
the dwarf galaxy scales (iBode et al. I (|200lh I. This effect 
reduces the expected abundance of dwarfs with respect to 
the predictions of CDM mod els, providing a better agr ee- 
rnent with observations fe.g.. iPapastergis et al. I (I2014ID . 

In this context, an attractive way to probe different DM 
cosmological scenarios is to look for faint dwarf galax¬ 
ies in the local universe and in particular in the Virgo 
cluster, which represents the densest nearby region. This 
large galaxy concentration allows us to probe the abun¬ 
dance of_dwarf_ga!axies_down to DM masses ~ 10 7 M 0 
(e.g.. ISabatini et al. (12005H 1 where pred ictions of WDM 
and CDM diverge ( Schneider et al. I(|2012l) h Nevertheless, 
the search for very faint dwarfs in Virgo is not straight¬ 


forward, as it implies pushing the detection toward very 
low limits in central s urface brightness, as first noted by 
ISandaee et al. I (|l985h . 

Impev et al. I (jl988H found a population of blue LSB 
dwarfs in the Virgo cluster with Mb < —11 not detected 
in previous surveys. They stressed the effects of this new 
LSB dwarf population on the faint-end of the Virgo galaxy 
luminosity functi on (LF hereafter). 

More recently iTrentham fc TulhH (l2002f) performed a 
survey in five different local fields of varying galaxy density 
including the center of the Virgo cluster. They obtained 
relatively deep Subaru images in a small area of the Virgo 
cluster surrounding the M86 bright galaxy. T hey made 
also u se of part of the ’VCC’ sample of iBinggeli et al. 1 
(|l985lf to derive the bright end of the LF. Their new Sub¬ 
aru dwarf sample allowed the evaluation of the Virgo LF 
down to M r ~ —12. Since the slope found at the faint end 
(—1.2) was significantly flatter than the CDM mass func¬ 
tion slope (—1.9), they suggested a high degree of suppres¬ 
sion of star formation in small DM halos. The suppressed 
luminosity of low-mass DM clumps would then explain the 
fl at slope of the luminosi ty function. 

ISabatini et al. I (|2003ll performed a wide survey (14 
deg 2 ) down to Mb ~ —10 with the Isaac Newton Telscope 
(2.5 m). Despite some success in discovering faint LSB 
dwarfs, they could not reach a high level of completeness 
and their derived LF appears flat in the magnitude interval 
—13 < Mb < —11. However, when combining their sam¬ 
ple with the brighter Virgo cluster sample, they suggest a 
somewhat steeper faint end LF than previously found. 

A recent investigation of the dwarf population us- 
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ing_CFHT imaging has been performed by iLieder et al. 1 
(1201^ , reaching a completeness limit of My ~ —13 in 3.75 
deg 2 . They found an average slope of the LF as steep as 
a = — 1.5 in the r ange —19 < M y < —13. 

Very recently, iDavies et al. I (120151) have provided a 
first analysis of the LSB dwarf population in a field of 
about 100 deg 2 f r om th e Next Generation Virgo Survey 
(iFerrarese et al. I (l2012fll. This survey is larger with re¬ 
spect to the one bv ISabatini et al.1 1 2003( 1 although not 
deeper, reaching completeness at g ~ 19 (which approx¬ 
imately corresponds to r ~ 18.5 for an average g — r ~ 
0.5). Their LF slope is consistent with previous values 
a ~ -1.35. 


To summarize, a population of faint LSB dwarfs in the 
Virgo cluster has emerged from imaging analysis, although 
the images were not deep enough to reach firm conclusions 
about the shape of the very faint-end of the Virgo LF. 

In this paper we show the results of a pilot LSB dwarf 
search in a deep Virgo field. Although smaller than pre¬ 
vious surveys, our observations are significantly deeper 
thanks to the use of th e Large Binocular Camera (LBC, 
iGiallongo et al. I (|2008lB of the La rge Binocular Tel e scope . 

The i mages were acquired by Lerchster et al. I (|2011h 
(see also iFassbender et al. I (|201l B in the context of the 
study of the high redshift [z ~ 1) galaxy cluster XMMU 
J1230.3+1339, that is located in the projected region of 
the Virgo cluster including NGC 4477. Although the field 
of view covers an area of only 0.16 deg 2 , it is particularly 
deep and provides a first estimate of the dwarf surface 
density at the faint-end of the galaxy LF in Virgo. 

We adopt = 0.7, flo = 0.3, and h = 0.7 in units of 
100 km/s/Mpc. AB magnitudes have been adopted except 
for Figure 5. 


2. THE LSB DWARF GALAXY SAMPLE 

The images were obtained in the U,B,r,i,z bands with 
exposure times ranging from 3.2 ks to 9 ks, and calibrated 
using zeropoi nts derived from Landol dt standards or Sloan 
DR7 catalog lAbazaiian et al. I (|2009h . The size of the field 
of view at full depth is ~ 57 6 arcmin 2 . Furth e r obs erva- 
tional details can be found in iLerchster et al. I (1201 lh . 

To detect faint LSB dwarfs and to confirm their associa¬ 
tion with the Virgo cluster we have used their distribution 
in th e uo — r plane, as previ o usly done by several au thors 
(e.g. IConselice et al. I (120021 ) . IRines fe Geller I (120081) and 
reference therein). We have performed the dete c tion i n 
the r band using SExtractor dBertin fc Arnoutn (|1996l )) 
setting a surface brightness limit of p ~ 27.2 at the 2er 
level and a minimum detection circular area of diameter 
2.5 arcsec. SExtractor has also been used to derive a first 
estimate of the central surface brightness. We plot the 
galaxies found down to r = 23 as a function of the central 
surface brightness in Figure 1. There are 11 galaxies that 
are well separated from the bulk of the galaxies in the cata¬ 
log (note that magnitudes and central surface brightnesses 
of the blue circles are more accurately derived from profile 
fitting, see below). A visual inspection has been executed 
to confirm the reality of such detections. All 11 LSB dwarf 
candidates are clearly extended sources with low surface 
brightness, and they are all detected in the full U,B,r,i,z 
multicolor dataset used. A further visual inspection over 
the whole image failed to detect additional candidates. 


To perform a more quantitative selection, we have plot¬ 
ted in the same Figure 1 a straig ht line represent i ng the 
LSB selection threshold used by IRines fe Geller I (120081) 
to separate Virgo galaxies from the background sample, 
extrapolated to our fainter magnitudes. 9 out of our 11 
candidates are above or very close to their selection thresh¬ 
old, and the remaining 2 are in any case well separated 
by the l ocu s of the backg round sample. We note that 
the IRines fe Geller I (|2008) threshold is a robust but not 
complete selection criterion, since there is a small overlap 
between the two populations, as shown by the spectro¬ 
scopic detection of a few Virgo galaxies below the selec¬ 
tion threshold (i.e. at brighter central surface brightness 
with respect to the adopted line for a given total magni¬ 
tude). As an independent check, we have also verified that 
dwarfs with similar characteristics are absolutely absent 
in our database of LBC fields even deeper than this (i.e., 
with exposure times ~ 3 — 10 hr and limiting magnitudes 
r ~ 26 — 27) obtai ned in different s ky re gions for a total 
area of ~ 0.6 deg 2 dBoutsia et al. I (120141) . Grazian et al. 
in preparation), lending further support to the association 
of the selected LSB dwarf candidates to the Virgo cluster. 
For these reasons, in the following we will consider all the 
II LSB dwarfs as bona fide Virgo members. 

The positions of all candidates in our image are shown in 
Figure 2 and their r-band morphologies are shown in Fig¬ 
ure 3. In the latter, their low central surface brightness 
with respect to the background sources is clearly evident. 
In some cases compact sources are overlapping with our 
LSB dwarf candidates, possibly altering the estimate of 
their total magnitude and morphological parameters. To 
obtain cl eaner measurements we have used the Galfit pack¬ 
age (jPeng et al. I (I2010I 1L a versatile software that allows 
radial profile modeling wit h differ ent functions. We used 
Sersic power-law profiles (jSersic I (|l968l) h without trun¬ 
cation, usually adopted to fit spiral as well as spheroidal 
morphologies. In the fitting procedure we leave the Sersic 
power-law index n, the scale radius r s , the axial ratio, 
the total magnitude (or the central surface brightness) as 
free parameters. The De Vaucouleurs and exponential disk 
profiles are included as cases with Sersic index n = 4 or 
1, respectively. The PSF profile was convolved with the 
intrinsic profiles for the fitting. Small, compact sources 
superimposed on our LSB dwarfs were simultaneously fit 
with the main LSB galaxy to remove their contribution. 
The most striking example is given by “K” galaxy which 
is contaminated by several compact sources, and which is 
shown in more detail in Figure 4. In this specific case 
differences with respect to the Sextractor parameters can 
reach values of the order of 0.4 in total magnitude. In 
general, however, differences between Galfit and Sextrac¬ 
tor total magnitudes or central surface brightnesses are no 
more than 0.15 mag or mag arcsec -2 , respectively. 

The resulting parameters, including magnitudes and 
central surface brightnesses, for the LSB dwarf candi¬ 
dates are shown in Table 1. For most sources the central 
surface brightness in the r band lies within the interval 
go ~ 24 — 26 and the Sersic scale radius in the interval 
r s = 3 — 10.7 arcsec or equivalently r s = 250 — 850 pc 
adopting an angular scale of 80 pc arcsec -1 , that corre¬ 
sponds to a distance of the Virgo cluster of 16.5 Mpc. 
Absolute magnitudes derived from the fit are in general 
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fainter than M r > —13 and reach values as faint as 
M r ~ —9.3 in the r band. Also, the Sersic index n is 
always n < 1 typical of high surface brightness spirals and 
dwarfs. In general the axial ratio is larger than 0.6, and 
for ~ 50% of the sample it is > 0.8. Formal statistical 
errors provided by Galfit are typically a few percent. 

To estimate galaxy colors in a similar way, without the 
contribution from overlapping compact sources, we have 
repeated the same procedure in our U,B,i,z images. We 
have again used Galfit, fixing the radial profile obtained in 
the r band image and leaving only the total magnitude free 
to change. In the estimate we assumed negligible color gra¬ 
dients along the galaxy profile. We have checked that the 
X 2 did not increase because of this assumption. The pro¬ 
cedure is often needed in the low signal-to-noise blue and 
UV images where only the central part of the galaxy profile 
emerges from the background noise. Results are shown in 
Table 1, where all the galaxies appear red with an average 
(B — r) ~ 0.8. To compare with the Sloan colors we fit the 
observed colors to theoretical spectral energy distributions 
(SED) using the procedure described below. We computed 
the Sloan g magnitude from the resulting best-fitting spec¬ 
trum, yielding an average (g — r) ~ 0.5. The latter color 
is somewhat bluer than that of the brightest Virgo galax¬ 
ies in our field (g — r ~ 0.9, 0.8 for NGC 4477 and 4479, 
respectively). The LSB colors are also blu er th an predic¬ 
tions from CDM models (g — r ~ 0.7, Wei nmann et al.l 
(j2011l lh which include specific prescriptions for star for¬ 
mation and feedback for simulated brighter Virgo dwarfs 
with M r < —15.5. As suggested by the same authors, the 
model predicts redder colors mainly because it could over¬ 
estimate the star formation quenching in Virgo dwarfs due 
to environmental effects (e.g., ram pressure stripping), re¬ 
quired to make their number density consistent with the 
observations. 

The LSB d warfs found in our field are similar to dwarfs 
fo und, e.g., bySabatini et al. I (l2003h . iLieder et al. I (|2012f) 
or iDavies et al. I ~ 2015f) but are considerably fainter. In 
partic ular we have only one dwarf in common with 
IDavies et al. I (l2015f ): in our data the ”K” dwarf (LSB- 
VCC025) appears more diffuse with a scale radius twice 
larger. This is a combined effect of the larger depth of our 
data and of our Galfit analysis, which shows that appro¬ 
priate photometry can be derived only after removing the 
faint background sources present inside the LSB profile 
(see Figure 4). For this reason, a direct com parison be¬ 
tween the magnitudes derived by us and by IDavies et al. I 
(1201511 must take into account the contribution by the faint 
interlopers, which amounts to Ar ~ —0.4. By assuming 
an average g — r ~ 0.4 (as derived from our table 1 with 
appropriate interpolation and filter co n versio n) we obtain 
the value g ~ 19.1 as in IDavies et al. I (1201511 . 

We have derived physical quantities for the 11 LSB 
dwarfs fitting the UBriz multicolor data to the syn¬ 
thetic SEDs based on a grid of models from the 
iBruzual &; Chariot! (j2003l f (BC) spectral synthesis code. 
The synthetic spectra are characterized by exponentially 
declining star formation histories of timescale r and by 
a set of ages, metallic i ties and dust ex t inctio ns (see, e.g., 
iGrazian et al. I (|2006lf , iFontana et al. I (1200611 for details). 
The redshift has been fixed to that of the Virgo cluster 
(peculiar velocities are irrelevant in this context). 


The output of this SED fitting procedure is shown in 
Table 1. All the galaxies are consistent with a stellar 
metallicity 0.2 solar and no dust. The resulting star for¬ 
mation rates (SFR) and stellar masses are indicative of old 
stellar populations with low SFR. In particular the ratio 
( age/r) ~ 3.6 is indicative of a passively declining stellar 
population, although this ratio is n ot as l arge a s in brighter 
ellipticals ((age/r) > 6 lFontana et al. I d2009j) ). 

Moreover, extrapolating the linear relation between star 
formation and stellar mass found for_ t he m ai n sequ ence 
galaxies in the local universe (see lElbaz et al. I (|2011l ff to¬ 
ward low masses, our dwarfs show average SFRs a factor 
10 lower for a given stellar mass with respect to the linear 
relation, consistent with a quenched star formation phase 
with typical specific SFR ~ 10 -11 yr -1 . 

3. THE VIRGO LUMINOSITY FUNCTION 

The main outcome of the present study is the large num¬ 
ber density of LSB dwarfs found in a single LBC pointing, 
corresponding to 11 LSB Virgo candidates (after remov¬ 
ing bright extended objects) in an effective field of 0.12 
deg 2 . The derived projected number density can be com¬ 
pared with corresponding proj ected densities of b r ighter 
Virgo galaxies. We use here the lTrentham fe Tullv I (l2002f ) 
sample, whose projected LF is computed per square mega- 
parsec, per magnitude and is normalized at a distance from 
the Virgo center (i.e. M87 position) of 200 kpc. We split 
our small sample into two bins (each 2 mags wide, one 
centered at M r = —12.5 and one at at M r = —10.5) that 
contain 4 and 7 candidates, respectively. 

Adopting for Virgo the angular scale of 80 pc arcsec -1 
we obtain a density of 200 and 350 objects Mpc -2 mag -1 
in our two bins at the position of our field, which is about 
1.3 deg or 340 kpc fro m M87. To rescale to 2 00 kpc 
for comparison with the iTrentham fc Tullv~l (I2002T) LF we 
have adopted their projected radial density distribution of 
brighter Virgo galaxies. We derive an increase by a factor 
~ 1.6 from 340 to 200 kpc. The resulting projected densi¬ 
ties of LSB dwarfs at 200 kpc scale to 320 and 560 Mpc -2 
mag -1 , respectively. 

These two average densities derived in our relatively 
small field ar e shown in Figure 5 togeth er with LF data 
adapted from ITrentham fe Tullv I (120021b Under the hy¬ 
pothesis that the small LSB dwarfs are distributed follow¬ 
ing the smooth gravitational potential well of the Virgo 
cluster we assume that the derived density is a fair repre¬ 
sentation of the average LSB dwarf density at ~ 300 kpc 
from the Virgo center. We explicitly neglect any possi¬ 
ble dwarf clustering or anti-clustering around the bright 
Virgo galaxies. Our sa mple shows an excess at M r > —17 
compared to the LF of ITrentham fe Tullv I (120021) but it 
is difficult to ascertain whether this is due to a re al devia¬ 
tion of the LF from a Schechter shape (ITrentham fc Tullv I 
(HU) or to observational biases which produce devia¬ 
tions from a featureless faint-end power law (see, e.g., 
ILieder et al. I (120121) h In practice the LF is de rived by 
connecting points from lTrentham fc Tullv I (120021) brighter 
than M r = —17, which are mainly based on the VCC sur¬ 
vey iBinggeli et al. I (|1985f) b with our two faint bins. The 
resulting slope is a ~ —1.4 and it is shown in Figure 5. 
We also show in the same figure a flatter —1.3 and steeper 
— 1.5 slope for comparison. Our value is consistent with 
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very rec ent estimates derive d from the overall Virgo re¬ 
gion bv iDavies et al. I (I2015T) a ~ —1.35, in the g band. 
However, their slope evaluation is limited to g ~ 18 which 
corresponds to r ~ 17.5 assuming g — r ~ 0.5. Although 
statistically limited by the small LBC area, our comple¬ 
mentary data allows us to probe the extension of the LF 
slope about 4 magnitudes fainter. 

4. THE IMPLICATIONS FOR DARK MATTER SCENARIOS 

The presence of a significant density of faint dwarfs in 
a rich cluster like Virgo can put interesting constraints on 
the nature of Dark Matter and on the physical mechanisms 
regulating galaxy star formation activity. In DM scenar¬ 
ios, galaxies form at the center of DM halos collapsed in 
over-dense regions. Following the progressive inclusion of 
previously collapsed halos into larger and larger structures, 
the central massive galaxies are then surrounded by a num¬ 
ber of satellites whose abundance depends on the nature 
of DM and on the baryonic physics driving the process of 
galaxy formation. 

CDM models tend to predict a large number density of 
small-mass ga laxies (M* < 10 lo L sign ifi cantly higher than 
is ob served dFontana et al. I (1200611 . iMarchesini et al. 1 
(2009)). In particular CDM models predict faint-end 
slopes of the galaxy LF that are much steeper than ob¬ 
served unless strong feedback from star formation activ¬ 
ity is assumed, particularly for the dwarf galaxy popu ¬ 
lation (ISomerv ille fc Primack I (llDDDh . I Cole et al~l ( 2000 1. 

(|2002tl l~ The feedback should be strong 


iMenci et al. 


enough to inhibit all star formation in many small CDM 
halos as recently suggested by ISawala et al. I (|2015t) . In 
this case, simple abundance matching between simulated 
DM halos and observed galaxies would not apply. 

Alternatively, the power of the DM spectrum at small 
scales can be reduced if density perturbations are damped 
below some characteristic scale by the action of the free 
streaming of WDM particles with therm al relic mass of the 
order of ~ 1 — 3 keV, see e.g.. ISommer-Larsen fc Dolgov I 
(j200ll ). This would decrease the predicted number den¬ 
sity of small DM halos, reducing the tension with the ob¬ 
served dwarf number density. Note that the suppression 
of the power spectrum with respect to the CDM case di¬ 
rectly depends on the particle mass m .x if the candidate 
is a thermal relic. For sterile neutrinos, the same sup¬ 
pression is obtained for a mass rn s teriie ~ 4.3 mx if the 
sterile neutrinos are prod uced f rom oscillations with active 
neut rinos in the Dodelson fc Widrow I (I1994T ) scenario (see 
IViel et al~ (12005)). For different production mechanisms 
(see Kusenko I (|2009lD the c onversion factor m^t erUe/mx 
ranges between 2 and 4 (see lDestri et al. I (I2013ID . 

Here we discuss the potential implications of our find¬ 
ings, and especially the constraints on the mass of the 
WDM particles. In particular, we compare the faint-end 
slope of the LF shown in Figure 5 with those derived from 
the sub-halo mass function in different DM models, includ¬ 
ing a proper treatment of the present uncertainty on the 
M/L ratio of dwarf galaxies. The mass dependence of the 
cumulative mass function can be described by a univer- 

1 OLDM 

sal form N{> M)/N to t oc [M/Mhost , independent of 

the mass of the host structure Mf lost , with an exponent 
determined by the assumed power spect rum. CDM simu¬ 
lations yield cxdm = —1-9 (see also lKlypin et al. 1 (|201llH . 


while smaller slopes are derived for WDM particles with 
mass mx ~ 1 keV (see, e.g., Knebe et al. 2005; Maccio & 
Fontanot 2010; Lovell et al. 2014). 

To derive the slope for different WDM models we need to 
explore the effects of different power spectra (correspond¬ 
ing to different DM particle masses) on the dynamical evo¬ 
lution of sub-halos. While full high-resolution N-body sim¬ 
ulations provide complete theoretical tools to perform such 
computations, they are also extremely time consuming and 
are generally applied to study specific DM models. In ad¬ 
dition, they are limited in resolution to sub-halo masses 
M/Mhost ^ 10 -4 . To effectively explore a wide range of 
DM (thermal relic) candidate masses (mx = 1, mx = 1-5, 
mx = 2 keV, and the CDM case mx 1 keV) with a res¬ 
olution extending down t o M/Mhost ~ 10~ 5,5 , we a dopt a 
semi-analytic model (see ISomerville fc Davdl (I2014T) for a 
review of N-body and semi-analytic models as theoretical 
tools to study galax y formation). 

This model (see IMenci et al. I (|2012f) . IMenci et al. I 
(2014J) for details) follows the merging histories of DM ha¬ 
los and sub-halos through a Monte Carlo simulation of 
the collapse and subsequent merging history of the peaks 
of the primordial density field. This lets us generate a 
synthyetic catalogue of model DM halos and of their past 
merging histories. Although the model can include a spe¬ 
cific treatment of the baryonic physics inside the halo’s 
evolution (e.g., gas cooling, star formation etc.), we prefer 
to derive constraints on the WDM particle mass assum¬ 
ing a general parametric power-law dependence between 
the galaxy luminosity and the total DM mass. The model 
includes the main processes affecting the dynamical evolu¬ 
tion of sub-haloes within larger DM haloes, i.e., dynamical 
friction, merging, and stripping . For a recent compar¬ 
ison of the dynamics of sub-halos in semi-analytic mod¬ 
els an d N-body si mulations see ISomerville fc Dave I (I2014T) 
an d lPullen et al. I (120141) . 

We first tested our semi-analytic modeling of the sub¬ 
halo population against previous works in specific test 
cases where detailed results are available. In the left 
panel of Figure 6 we have compared the sub-halo mass 
functio n obta i ned f rom our model with those obtained by 
I Pullen et al. I (l2014f) for a CDM spectrum and for a WDM 
spectrum corresponding to a thermal relic particle mass 
m,x = 1.5 keV. Those authors adopt a semi-analytic model 
calibrated and tested through the comparison with ultra 
high-resolution N-body simulations. The average slopes 
of the sub-halo mass functions we derive match the exist¬ 
ing results not only for the CDM case (black curve, black 
points) but also for the WDM case (red dotted curve, red 
points) with mx = 1.5 keV. 

As we increase the relic WDM particle masses in our 
model, we obtain increasing slopes for the DM sub-halo 
mass distribution. In particular we obtained aoM = —1.2 
for mx = 1 keV, odm = —1-4 for mx = 1.5 keV, 
oldm = —1-6 for mx = 2 keV, and <xdm = —2 for the 
CDM case. 

To convert the sub-halo mass distribution N(M) into 
a LF (i.e. a N(L)) we parametrize the mass-to -light 
ratio at low masses as M/L oc M x ~$ (jMoster et al. I 
(120131) ). so that the slope of the LF can be written as 
a = (1 + ex dm ~ P)//3 assuming proportionality between 
the galaxy luminosity and stellar mass. This turns out 
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to be a good approximation for a limited range of (faint) 
magnitudes and for fairly old stellar populations. Both 
assumptions are satisfied here according to our SED fit¬ 
ting. Since at present substantial uncertainty affects the 
measured M/L ratio for dwarf galaxies, our strategy is 
based on bracketing the present uncertainty by explor¬ 
ing values of 0 in a wide range from 1 to 3. A value 
of 0 > 1 is der i ved fr om abu ndance matching techni ques 
(|Moster et al. I (l2013t) : see also lBehroozi et all (12013IB as¬ 
suming a standard cuspy (e.g. Navarro, Frenk & White) 
inner density CDM profiles for the observed satellites in 
the Local Group. In this case only a fraction of CDM ha¬ 
los can host galaxies. More complex and shallower density 
profiles assuming a dynamical interaction between baryons 
and DM halos are consis t ent w ith a shallower slope 0 ~ 1 
(see iBrook fe Di Cintiol (1201511 and ISawala et al. I (1201511 
for both cases). 

Assuming that all halos host a galaxy, we can then 
derive the slope a of the LF obtained in different 
CDM/WDM scenarios as a function of the assumed value 
for 0, and compare it with our measured value a = 
— 1.4 ± 0.05 derived from our dwarf sample. The output 
of this exercise is shown in the right panel of Figure 6. 
We confirm that the CDM model needs a relatively large 
value of 2 < 0 < 2.5 to be compatible with the data. 
In the case of WDM models, we are able to set a lower 
limit of nix > 1-5 keV to be consistent with the mea¬ 
sured slope a , adopting lower values of 1 < 0 < 2. On 
the contrary, WDM models with mx < 1-5 keV are not 
able to provide the observed faint-end LF slope if 0 = 1 
represents a reasonable lower limit for the M/L relation. 
Of course, our low limit on the WDM particle mass could 
even increase if we allow for a complete inhibition of star 
formation in a fraction of the low-mass sub-halo popula¬ 
tion due to reionization and backg round UV r adiation, as 
suggested for instance bv ISawala et al. 1 (|2015f l. 

This limit compares well with other existing limits. 
Lower limits mx > 1 keV have been derived from galaxy 
counts in the Hubble Ultra Deep Field (ISchultz et al. 1 
(12014H L while limits mx > 1-5 keV have been derived from 
the abundance of ultra-faint Milky Way sat elli tes mea¬ 
sured in the Sloan Digital Sky Survey (see iHoriuchi et al~l 
(2014]) and references therein). Much tig hter constraints, 
m x > 3.3 keV (e.g. IViel et al~l (12013IB . have been de¬ 
rived by comparing the observed Lyman-a forest of high- 
resolution (z > 4) quasar spectra with hydrodynamical 
N-body simulations. Various uncertai nties may sti ll af- 
fect these constraints ( see di scussions in lAbazaiian et al~l 
(|201lh . ISchultz et al. I (l2014l) h The comparison of subha¬ 
los to Milky Way dwarfs assumes a factor « 4 correction 
for the number of dwarfs being missed by current surveys, 
and lower correction factors would appreciably weaken the 
constraints. Lyman-a absorption is also a challenging tool 
and requires disentangling the effects of pressure support 
and thermal broadening from those caused by DM spec¬ 
trum, as well as assumptions on the thermal history of the 
intergalactic medium and of the ionizing background. 

5. SUMMARY 

We have used deep LBC exposures obtained in five 
bands (UBriz) over a 576 arcrnin 2 pointing located within 
the Virgo field to find and study faint LSB dwarfs asso¬ 


ciated with the Virgo cluster at an unprecedented depth 
of M r ~ —9. Despite the relatively small field of view, 
we detect in the r band as many as 11 LSB dwarfs in the 
magnitude interval —13 < M r < —9, which we identify as 
Virgo members on the basis of their position in the go — r 
plane. 

Due to the significant depth of our images, their mor¬ 
phology appears to be contaminated by smaller, or com¬ 
pact sources, most likely in the background, which would 
affect the estimate of their total flux and morphological 
parameters. To remove them from our analysis, we have 
used Galfit to separately model the LSB dwarfs and the 
other sources, in order to derive contamination-free mor¬ 
phological parameters and magnitudes. This procedure 
has been applied to all images in five bands. We find that 
our sources have scale radii 250 < r s < 850 pc and an 
average red color U — r ~ 2.1. Based on this evidence, we 
claim that these sources represent the faintest extremes 
of the population of low luminosity LSB dwarfs that has 
been recently detected in wider s ur veys of the Virgo cluste r 
(e.g. Trentham fc Tullv (l2002f). IS abatini et al. I (|2003l) . 
iLieder et al. I ( 201211 and Davies et al. I ( 201511 1. 

The average ratio between the age of the last episode of 
star formation and the exponentially declining star forma¬ 
tion timescale (age/r) =3.6 indicates a passively evolving 
phase with a specific SFR of ~ 10 -11 yr -1 . They are older 
(~ 7 Gyr) than the typical cluster formation age and ap¬ 
pear regular in morphology. 

We have converted our detections into average pro¬ 
jected number densities, normalized at 200 kpc from the 
Virgo center, computed in two magnitude bins centered at 
(M r ) = —12.5 and (M r ) = —10.5. We find a number den¬ 
sity of 320 and 560 sources Mpc -2 mag -1 , respectively. 
These results allow us to extend the estimate of the Virgo 
LF down to very faint limits. Merging our data with the 
brighter points mainly based on the VCC survey, we find 
an average faint-end slope of about a czl —1.4 down to 
M r -9.5. 

It is well known that the relatively shallow slope of the 
LF sets challenging constraints on Dark Matter scenarios. 
Strong baryonic feedback processes must be invoked to rec¬ 
oncile the intriniscally steep CDM slopes with the observed 
flatter LF. Instead, we explore the alternative options, i.e., 
those provided by WDM scenarios with particles of ~keV 
mass scale._ 

Using our semi-analytic model (|Menci et al. Il2012i f we 
show that the slope ~ —1.4 we derive for the LF can be 
translated into an interesting low limit on the WDM parti¬ 
cle mass >1.5 keV, regardless of the exact relation present 
at low masses between the M/L ratio and the halo mass 
M (i.e. regardless of the details of the baryonic feedback 
processes). This limit compares well with previous limits 
derived from gal axy co unts in the Hubble Ultra Deep Field 
(jSchultz et al. I (I2014I 1L from the abundance of ultra-faint 
Milky Wav sate l lites i n the Sloan Digital Sky Survey (see 
IHoriuchi et al. I(l2014fll or f rom the observed Lyman-a for¬ 


est (e.g., IViel et al. 120131 U. 


It is intriguing to see that this small but non-negligible 
sample has been detected in a single pointing taken with 
a ~ 24' x 24' imager at an 8m telescope (LBC at LBT 
in this case). While it is still limited by the small area 
surveyed so far and by the lack of spectroscopic follow up, 
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the extension of the same analysis to a larger Virgo region 
using deep (2-4h of exposure time per band) multicolor 
images at 8m class telescopes is clearly needed to reduce 
the current uncertainties. 

We thank the referee for useful comments which have 
significantly improved our paper. We also thank M. Dick¬ 
inson for critical reading. Observations have been carried 
out using the Large Binocular Telescope at Mt. Graham, 
AZ. The LBT is an international collaboration among in¬ 


stitutions in the United States, Italy, and Germany. LBT 
Corporation partners are The University of Arizona on be¬ 
half of the Arizona university system; Istituto Nazionale di 
Astrofisica, Italy; LBT Beteiligungsgesellschaft, Germany, 
representing the Max-Planck Society, the Astrophysical 
Institute Potsdam, and Heidelberg University; The Ohio 
State University; and The Research Corporation, on behalf 
of The University of Notre Dame, University of Minnesota, 
and University of Virginia. 
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Table 1 

LSB CATALOG IN THE VlRGO-XMMUJl230 FIELD 


ID 

A 

B 

C 

D 

E 

F 

G 

H 

J 

K 

L 

RA (187) 

0.3919 

0.6635 

0.5879 

0.5483 

0.3814 

0.3863 

0.6600 

0.6224 

0.4138 

0.7070 

0.7244 

DEC (+13) 

0.7704 

0.7391 

0.7058 

0.6914 

0.6619 

0.6222 

0.6196 

0.5565 

0.5075 

0.5053 

0.4939 

r 

21.8 

20.4 

21.2 

18.3 

18.8 

21.9 

18.1 

21.2 

21.2 

19.1 

21.9 

M0 

24.9 

24.5 

25.8 

24.0 

23.5 

24.5 

24.3 

25.2 

25.5 

25.8 

26.1 

n a 

0.9 

0.7 

0.7 

0.7 

0.7 

0.9 

0.6 

0.5 

0.7 

0.4 

0.3 

M r 

-9.3 

-10.7 

-9.9 

-12.8 

-12.4 

-9.3 

-13.1 

-9.9 

-9.9 

-12.0 

-9.3 

r s (arcsec) c 

3.1 

4.5 

6.1 

8.0 

5.7 

2.7 

10.4 

3.0 

4.2 

10.7 

3.0 

r s (pc) c 

248 

360 

488 

640 

456 

216 

832 

240 

336 

856 

240 

axis ratio 

0.7 

0.7 

0.7 

0.9 

0.7 

0.6 

0.7 

0.9 

0.8 

0.9 

0.9 

r — i 

0.1 

0.1 

0.2 

0.2 

0.1 

0.1 

0.1 

0.3 

0.2 

0.0 

0.0 

i — z 

0.2 

0.2 

0.0 

0.2 

0.2 

0.0 

0.2 

0.2 

0.0 

-0.2 

0.1 

B-r 

0.8 

0.9 

0.7 

1.0 

0.9 

0.9 

0.9 

0.8 

0.9 

0.8 

0.7 

U - B 

0.8 

1.0 

1.4 

1.2 

1.0 

1.2 

1.1 

1.2 

1.5 

2.2 

0.8 

SFR d 

2 

5 

2 

17 

22 

1 

30 

2 

0.1 

37 

4 

M* e 

0.5 

2 

1 

20 

10 

0.5 

22 

1 

1 

4 

0.3 

ages-f 

7 

7 

7 

9 

7 

7 

8 

8 

8 

6 

6 

t 9 (Gyr) 

3 

2 

2 

2 

2 

2 

2 

2 

1 

3 

5 


a Sersic index _ 

b Absolute magnitudes computed adopting an average distance modulus for Virgo AM = 31.1 dMei et al. I <120071) ). 
c Scale radius from the Sersic profile fitting. An angular scale of 80 pc arcsec -1 has been adopted. 

d Star formation rate in units of 10 —5 Mq yr —1 . Resulting best fit Bruzual & Chariot models have no dust and 0.2 solar metallicity. SFR are 
uncertain by a factor 1.2 on average. 

e Stellar mass in units of 1O 6 M0. Average uncertainties are of the order of 40%. 
f Ages in units of Gyr. The 68% average probability distribution ranges from 6 to 12 Gyr. 

9 Exponential timescale of declining star formation. The 68% average probability distribution ranges from 2 to 8 Gyr. 












Fig. 1.— Surface brightness - magnitude relation for the galaxies in the field. Faint sources located in the extended halos of bright galaxies 
and blended galaxies have been removed from the plot since they have biased photometry. LSB dwarf candidat es are shown as b l ue ci rcles. 
The continuous line represents an extrapolation to fainter magnitudes of the LSB selection threshold adopted bv IRines fc Geller ( f 2008h 
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Fig. 2.— Position of the selected LSB dwarfs in the Virgo-xmmujl230 LBC field 
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Fig. 3. — Selected LSB dwarfs in the Virgo-xmmujl230 field; the box size of each image is ~ 57 arcsec. The sequence from the top-left to 
the bottom follows the list in Table 1. 
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Fig. 4.— The ”K” dwarf galaxy. From left to right: observed profile; model profile from Galfit (the brightest background galaxies have 
been fitted separately); residuals. 
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Fig. 5.— Virgo projected LF normalized at 200 kpc. Filled squares indicat e data from the present sa mple after conversion from AB to the 
Vega magnitude system by R(V ega) = r(AB) — 0.2. Empty squares are from I Trent ham Sz Tullvl d2002D in the Vega system. The continuous 
curve represents a Schechter shape with slope a ~ —1.4 and M* ~ —22.5. Two faint slopes a ~ —1.3, —1.5 are also shown for comparison 
(dashed and dotted lines, respectively). 
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Fig. 6. — DM model predictions. Left: Cumulative mass fun ction vs satell i te ove r main halo mass ratio predicted by our Monte Carlo 
model (black and red curves) compared with results derived by I Pullen et al.H d2014T) . Black curve and points are for the CDM spectrum 
and red curve and points for WDM spectrum with rax = 1.5 keV. Right: luminosity (i.e. stellar mass) function slope a vs (3, the slope of 
the M/L oc M 1— P power-law relation. The observed value a = —1.4 =b 0.05 is shown as a dashed region. Different curves derived from our 
semi-analytic model are shown for different WDM particle masses and CDM (1,1.5,2 keV and CDM starting from top to bottom). For (3 > 1 
then rax > 1.5 keV if a ~ —1.4. 







